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Abstract 

It is shown that under essentially all conditions, the non-linear classical equa- 
tions governing gravitation and matter in cosmology have a solution in which 
far outside the horizon in a suitable gauge the reduced spatial metric (the 
spatial metric divided by the square of the Robertson- Walker scale factor 
o) is time-independent, though with an arbitrary dependence on co-moving 
coordinates, and all perturbations to the other metric components and to all 
matter variables vanish, to leading order in 1/a. The corrections are of order 
1/a 2 , and are explicitly given for the reduced metric in a multifield model 
with a general potential. Further, this is the solution that describes the 
metric and matter produced by single- field inflation. These results justify 
the use of observed non-Gaussian correlations (or their absence) as a test 
of theories of single-field inflation, despite our ignorance of the constituents 
of the universe while fluctuations are outside the horizon after inflation, as 
long as graphs with loops can be neglected. 
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I. Introduction 



Non-Gaussian cosmological correlations are attracting increasing interest 
as an observational test of detailed theories of inflation [1]. But there is a 
problem in calculations of observable non-Gaussian correlations. Given any 
specific Lagrangian for the scalar fields that play a role in inflation, we know- 
in principle how to calculate the correlation functions for these fields and 
gravitation up to the end of inflation. And given any set of correlation 
functions for gravitational and matter and radiation perturbations at some 
time in the relatively recent era when the temperature is well below the 
QCD scale, we know enough about the contents of the universe to calculate 
the subsequent evolution of these correlations. But in the intervening era 
the universe went through a sequence of transformations about which we 
know almost nothing, including reheating, lepton and baryon synthesis, and 
cold matter decoupling. So how can we use assumptions about inflation to 
calculate observable correlations? 

The only thing that gives us any hope in understanding cosmological 
correlations is that the wavelengths at which these correlations are observed 
were outside the horizon during the whole period from a time during in- 
flation until a relatively recent time when the contents of the universe are 
reasonably well understood. But in order to take advantage of this fact, 
we need to identify a set of variables whose correlation functions are time- 
independent for wavelengths outside the horizon. In the linear approxima- 
tion, it is known that quantum fluctuations in single-field inflation produce 
adiabatic fluctuations, in which the curvature perturbation ( as well as the 
amplitude of gravitational waves become constant outside the horizon[2]. 
This is enough to show that Gaussian correlations of these quantities re- 
main time-independent after inflation, as long as the wavelength is outside 
the horizon, at least when quantum effects are neglected. But for non- 
Gaussian correlations we must work with the full non-linear field equations. 

It is known that the classical non-linear field equations for single-field 
inflation have an " adiabatic" solution for which ( and the gravitational wave 
amplitude become time-independent at late times during inflation[3], and it 
has further been shown[4] that these quantities become time-independent 
outside the horizon both during and after single-field inflation, but this 
is only part of the story. To provide initial conditions for a calculation of 
fluctuations through horizon re-entry and until the present, we need to know 
not only £ and the other metric components but also the matter (including 
radiation) perturbations before horizon re-entry. It is sometimes taken as 
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part of the definition of the adiabatic solution that (in a suitable gauge) 
these matter perturbations vanish, but it needs to be shown that such a 
solution exists, and that in some circumstances the universe is described by 
this solution. 

Section II of this paper shows that whatever the constituents of the uni- 
verse and the classical equations governing them may be, these equations 
have a solution for which in a suitable gauge, as long as all relevant wave- 
lengths are sufficiently far outside the horizon, all components of the reduced 
metric gij = gij /a 2 become time- independent functions of position; goo be- 
comes — 1; <7io vanishes; and all matter densities, pressures, and velocities 
become equal to their unperturbed values; in all cases with corrections of 
order (k/aH) 2 . (As usual, a(t) is the Robertson-Walker scale factor and 
H(t) = a(t)/a(t), while k is the largest relevant wave number.) The ar- 
gument for these results is based on considerations of broken symmetry, 
similar to those used to derive the form of the chiral Lagrangian for soft 
pions[5]. It relies only on the general covariance of the underlying equa- 
tions, and the usual assumption that these equations have a solution of the 
Robertson-Walker form. This argument may be regarded as a substitute for 
a "separate universe" assumption [6], but it gives more detailed information, 
and some may find it more convincing. 

In Section III we verify these results in a fairly general model of matter 
fields, in a special gauge that allows the calculation of metric perturbations 
outside the horizon without at the same time having to solve the equations 
for matter perturbations. In particular we confirm that in this model there is 
always an adiabatic solution in which the corrections to the leading terms for 
both the metric and the matter variables are of order (k/aH) 2 . This much is 
already known for the metric components [4], but here we also obtain explicit 
results** for the form of the 0(l/a 2 ) and 0(l/a 3 ) terms in <7ij[7]. In general 
this is only one of many possible solutions; whether or not it is the solution 
that describes the real world depends on the details of inflation. In Section 
IV we show that single field inflation in this model leads to matter and 
metric fields that are described by the adiabatic solution after inflation, as 
long as all wavelengths are outside the horizon. 

There is a problem with all such classical field calculations. Even if we 
could show that under all circumstances the full non-linear classical field 

"Added note: In a paper now in preparation, I show that these explicit results are not 
limited to scalar field theories, but are quite general for theories with anisotropic inertia 
and vorticity that vanish to order aT 2 . 
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equations have a solution for which ( and gravitational wave amplitudes be- 
come constant outside the horizon, and that the initial conditions provided 
by single-field inflation (or a state of thermal equilibrium after inflation) pro- 
duce perturbations that are described by this solution, we still would not 
know that Heisenberg picture quantum operators for ( and gravitational 
wave amplitudes become constant at late times, and so we would not know 
that the correlation functions become constant outside the horizon. "Out- 
side the horizon" means that for all relevant co-moving wave numbers k, 
we have k/a <C H. For quantum operators there can be no clear meaning 
to this, because whatever the wave numbers at which the correlation func- 
tions are observed, quantum fluctuations can carry arbitrarily large virtual 
wave numbers. That is, there is no limit on the wave numbers circulating in 
the loops in general graphs, however small we make the wave numbers for 
the external lines of these graphs. For inflation with a single inflaton field, 
the loop contributions to the correlation function of ( are integrals over 
virtual wave numbers p with integrands that are time-independent when 
the virtual as well as the external wavelengths are outside the horizon, but 
virtual wavelengths can not be constrained to remain outside the horizon, 
because the integrals over virtual wave numbers are ultraviolet-divergent. 
(For examples, see [8].) True, we can assume that the ultraviolet divergence 
is canceled by counterterms arising from \J — Detg W l 'R^ v and \J — Detg R 2 
terms in the Lagrangian density [9], but if this cancelation results in an ef- 
fective cut-off at p/a of order H, then the correlation functions will involve 
powers of In a [8]. Detailed calculations [10] confirm the presence of such 
time-dependent corrections in loop contributions to correlation functions. 

Fortunately, in many theories the tree graphs make much larger contri- 
butions to the correlation functions than loop graphs. For a tree graph the 
wave number associated with any internal line is just a sum of wave num- 
bers of several external lines, so all internal wavelengths can be assumed 
to be outside the horizon if the external wavelengths are. Thus at least in 
some theories, one can treat the non-linear field equations as if all relevant 
wavelengths were outside the horizon, and hope that quantum effects do not 
introduce large corrections. Alternatively, one can limit oneself to the tree 
approximation from the beginning, and assume that tree graphs give a good 
approximation to the correlation functions. 

A recent paper [11] showed how to calculate the sum of tree graphs for 
the generating function for general correlation functions by solving the clas- 
sical equations of motion subject to certain constraints that depend on the 
current appearing in the generating function. This is reviewed here in an ap- 
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pendix, using a simplified notation and adding some necessary comments. 
In order to conclude in this formalism that correlation functions become 
time-independent outside the horizon, it is not enough to show that the so- 
lution of the non-linear classical field equations becomes time-independent 
at late times during inflation. As reviewed here in the appendix, one must 
also show that the effect of the constraints that are imposed at the time at 
which the correlation functions are measured becomes independent of this 
time when all wavelengths are outside the horizon, and also that a certain 
integral converges. In Section V we show that these conditions are all sat- 
isfied during and after single-field inflation for the quantity = gij/a 2 . 
(The same argument applies to any function of g~ij, such as the quantity 
£ = In ^/Detg studied in [3].) Thus in order for parametric amplification 
during reheating[12] to produce significant changes in the correlation func- 
tions, such effects would have to amplify perturbations by a factor of order 
e 120 to e 140 . 

In summary, these results provide a practical program for calculating 
observable correlation functions from theories of single-field inflation. 

(i) First calculate the correlation functions of gij = gij/a 2 (or any functions 
of cjij) sufficiently late after horizon exit during inflation so that they are 
time-independent, using a definition of the time coordinate for which the 
inflaton field is unperturbed. (This would presumably be done by direct 
calculation of tree graphs, as already done in [3] for the bispectrum, rather 
than by using the methods of Section V, which are intended only to provide a 
proof of the time-independence outside the horizon of the sum of tree graphs 
for any correlation function of gij.) If we like we can separate a curvature 
perturbation ( by following [3] and writing 

9ij = e 2C [e%- , Tr 7 = 0, 
or alternatively by writing 

9ij = S ij + 2 C Sij , Tr 7 = . 

These different definitions of course give different non-Gaussian correlation 
functions for (, but with either definition the correlation functions are con- 
stant outside the horizon. 

(ii) At a time which is sufficiently early so that the all wavelengths are still 
outside the horizon, but late enough so that the contents of the universe are 
well understood, take the correlation functions of or of functions of gij 
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to be given by the results of (i), and take all correlation functions involving 
<7i0 and/or goo + 1 and/or matter or radiation perturbations to vanish. 

(hi) Use the results of (ii) as initial conditions for calculation of the subse- 
quent evolution of the correlation functions for gravitational and matter and 
radiation perturbations when the wavelengths re-enter the horizon. This can 
be done by using the classical field equations to derive coupled differential 
equations for the various correlation functions, but such calculations are 
outside the scope of this paper. 

The above program appears to be more or less what is done in recent work 
on non-Gaussian correlations [13]. The aim of this paper is to clarify the 
justification for these calculations. 



II. The General Adiabatic Solution 



We assume, as usual in cosmology, that whatever the dynamical equa- 
tions governing the metric and matter (including radiation) variables may 
be, these equations have a solution in which the metric takes the Robertson- 
Walker form, with goo = —1, goi = 0, and gij = gij/a? = 5ij, and in which all 
matter variables take their unperturbed form; that is, all densities and pres- 
sures and scalar fields are functions only of time, and all velocities and other 
3-vectors vanish. This section will give a very general argument that, what- 
ever the constituents of the universe and the generally covariant equations 
governing them and the metric may be, for a suitable choice of spacetime 
coordinates, these equations always also have a family of solutions that we 
will call adiabatic, which for large a(t) have the following properties: 

1. The metric for any of these solutions has components with 

S0o(x,t) = -l + 0(a- 2 (t)) , 5i0 (x,i) = o(a- 2 (t)) , 
gij {x,t) = a 2 {t) [g tj {x) + 0(a- 2 (t))] , (1) 

where Gij{x) is an arbitrary function only of the spatial coordinates. 
(Different choices of this function characterize the different members 
of this family of solutions.) 

2. Whether or not the energy and momentum of any particular con- 
stituent of the universe is separately conserved, its energy-momentum 
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tensor has the form 



T 00 (x,t) =p(t) + 0(a- 2 (t)) , T l0 (x,t) = 0(a- 2 (t)) 
Ty(x,t) = a 2 (t) [g tj (x)p(t) + o(a~ 2 (i)) . 



(2) 



(Here and below, a bar over any quantity indicates its unperturbed 



3. Any four-scalar s(x, t), such as a temperatures, number density, or 
scalar field, has the form 



We are not assuming a de Sitter expansion, but in counting powers of I /a, 
we shall take H and its time derivatives to be of zeroth order in a, so that 
quantities like a(t) and / a(t) dt are counted as being of first order in a. It 
should be understood that since the scale of a has a physical significance 
only when a multiplies a co-moving coordinate, it follows that when we 
calculate correlation functions with a typical co-moving wave number k, 
a factor a -1 will always be accompanied with a factor k. Since k/a has 
the same dimensions as H = a/ a, we can anticipate that the dimensionless 
parameter that characterizes the smallness of a term of order a~ n is {k/aH) n . 
Thus this theorem gives good approximations to the adiabatic solutions both 
after horizon exit during inflation, when k/aH is decreasing, and before 
horizon re-entry after inflation, when k/aH is increasing, as long as k/aH 
is sufficiently small. 

Of course, these solutions are in general far from unique, and the state- 
ment that these adiabatic solutions exist does not tell us that one of these 
solutions actually describes the metric and matter of the universe. As we 
will see in section IV, if we start with single-field inflation then the universe 
will thereafter be described by an adiabatic solution. Also, even when the 
universe is described by an adiabatic solution, we need a detailed model of 
inflation to calculate the function in Eq. (1). 

To prove the existence of the adiabatic solutions, we will make use of an 
argument based on the broken symmetry of general covariance. As already 
mentioned, we are assuming that the dynamical equations have a solution 
in which the metric takes the Robertson-Walker form, and in which all mat- 
ter variables take their unperturbed form, with pressures and densities only 
functions of time, and vanishing co-moving velocities. Now, whatever they 



value.) 




(3) 
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are, the dynamical equations will be invariant under all coordinate transfor- 
mations, but this solution is not. In particular, if we subject the space coor- 
dinates to a matrix transformation x l — > x H = A l jX^ , with A l j an arbitrary- 
constant real matrix, then we get another exact solution, with goo = — 1) 
gio = 0, but now with = gij/a 2 equal to the arbitrary constant positive 
real matrix (A T A)ij. The energy-momentum tensor of any constituent of 
the universe (whether or not separately conserved) will in the new coordi- 
nate system still have the perfect fluid form, T^ u = g^ u p+u ll u u (p + p), with 
the same density p(t), pressure p(t), and velocity Ui = 0, uq = — 1, but now 
with the new metric. 

Instead of this exact solution, now consider what we will call a "trial 
configuration" in which goo = — 1) 9i0 = and all densities, pressures, and 
velocities are unperturbed, but with an arbitrary time-independent pos- 
itive matrix function Gij{x) of the co-moving space coordinates x 1 , not nec- 
essarily close to 5ij . This trial configuration is of course not a solution of the 
field equations, but since it would be a solution if gij were constant, it fails 
to be a solution only because there are terms in the field equations in which 
space derivatives act on g^. (Up to this point, this is just like the argument 
used to derive the effective chiral Lagrangian for soft pions[5].) The spatial 
derivatives of Qij(x.) thus act as forcing terms, that drive the actual solution 
away from the trial configuration. That is, making the tentative assump- 
tion that the differences between metric or matter variables and their values 
in the trial configuration are small perturbations when a(t) is sufficiently 
large, these perturbations satisfy a set of coupled inhomogeneous linear dif- 
ferential equations, with left-hand sides that are linear combinations of time 
derivatives of these perturbations, and right-hand sides that involve spatial 
derivatives of Gij (x) . We will see concrete examples of such equations in the 
next section. 

Now, as a special case of general covariance, the field equations must be 
invariant under the substitution x % — > \x % (with A an arbitrary constant) if 
we also subject other quantities to appropriate transformations: 3-tensors 
such as g^ and Tjj transform as g^ — > \~ 2 gij and — > A~ 2 Tjj, while 3- 
vectors such as gio and Tjo transform as g^ — > X~ 1 g-io and T^o — > A -1 Tjo- 
(Here T^ u may be the energy-momentum tensor of any one constituent of 
the universe, even if not separately conserved.) It is convenient to express 
this as invariance under a scale transformation: 

x i Ax* , a(t) -► A _1 a(t) , (4) 
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that leaves invariant various reduced quantities 



9ij{x,t) = gij(x,t)/a 2 (t) , 
f i3 -(x,t) = T tl (x,t)/a 2 (t) , 



ffio(x,t) = g i0 (x,t)/a(t) , 
fiofct) = T i0 (x,t)/a(t) , 



as well as all 3-scalars such as temperature, densities, scalar fields, and also 
500 an d Tqq. The forcing term for any scale- invariant perturbation must be 
scale-invariant, and three-dimensional coordinate invariance requires it to 
have the same transformation under purely spatial coordinate transforma- 
tions as the perturbation, so the perturbation of any scale-invariant quan- 
tity away from its value in the trial configuration will be proportional to as 
many powers of 1/a as appear in the scale-invariant quantity formed from 
1/a and derivatives of Qij that has the same transformation property under 
three-dimensional coordinate transformations as the perturbation in ques- 
tion. For perturbations to the scale-invariant quantities gij/a 2 or Tjj/o 2 , 
the scale-invariant forcing terms with the minimum number of factors of 
1/a are proportional to the 3-tensors TZij/a 2 or QijlZ/a 2 , where IZij is the 3- 
dimensional Ricci tensor for the 3- metric Qij, and 1Z = Q kl 1Zki, with Q^ the 
reciprocal otQij. Likewise, for perturbations to the scale- invariant quantities 
gio/a or Tjo/a, the scale- invariant forcing terms with the minimum number 
of factors of 1/a are proportional to the scale-invariant 3- vectors dilZ/a 3 , 
and for perturbations to the scale- invariant quantities goo or Too, the forcing 
terms with the minimum number of factors of 1/a are proportional to the 
scale-invariant 3-scalar IZ/a 2 , and likewise for perturbations to any other 
scale-invariant 3-scalar. Thus the difference between the values of the quan- 
tities gij/a 2 , Tij/a 2 , g^, T^o, 500 , Tqq, and 3-scalars like temperatures or 
scalar fields and the values of the corresponding quantities in the trial con- 
figuration are all of order 1/a 2 , as was to be proved. In particular, all these 
perturbations are small for sufficiently large a{t), as tentatively assumed in 
proving the existence of these solutions. 

The general solution for the perturbations to the trial configuration con- 
sists of a sum of the solution of the inhomogeneous differential equations, 
with derivatives of Qij (it) as forcing terms, plus solutions of the correspond- 
ing homogeneous equations. In a completely general theory of inflation the 
solutions of the homogeneous equation could have any magnitude. However, 
we will see in Section IV that in single field inflation they are also of order 
a- 2 . 

The form (1), (2), for the adiabatic solutions is not valid for all choices 
of spacetime coordinates, but it is easy to impose gauge-fixing conditions 
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on the coordinates that are consistent with this form. We can choose the 
time-coordinate so that any one three-scalar, such as a scalar field or the 
temperature, is unperturbed. (A generalized version of this choice of gauge 
is adopted in Section III.) To choose the space coordinates, we note that 
under a time-dependent transformation x % — ► x ll (x, t) that leaves the time 
invariant, the metric component g t0 undergoes the transformation 

jo _^ ho _ ® x ' % jo , ^ x ' % „oo 

We can evidently choose the time-dependence of x /l (x, t) so that g' l ° = 0, 
by solving the differential equation 

dx* = _dx* , 00 
dt dxi y /y 

for any arbitrary choice of x H (x, to) at an initial time to- I n this case, also 
g' i0 = 0. Though not unique, this choice of space and time coordinates is 
clearly consistent with (1) — (3). It still leaves us free to make purely spatial 
time-independent coordinate transformations, a freedom we have used in 
the arguments above. 

These adiabatic solutions to the non-linear field equations far outside 
the horizon may look unfamiliar to readers who are familiar with the form 
of the adiabatic solution for scalar modes in the linear approximation in 
Newtonian gauge, for which g-iQ = and g-ij oc 5ij. In [2] it is shown in the 
linear approximation that in the adiabatic mode in Newtonian gauge there 
are perturbations to the matter fields Xn{x) that do not vanish for large 
a(t), and a perturbation to g^ that does not become time-independent in 
this limit: 



<5 Xn (x, t) = - C(X) ,X f) / a(f ) dt' , 5~ 9ij = 2<%C(x) 1 - f a(t>) dt' 

a(t) Jt V a(t) Jt 

with C(x) an infinitesimal function only of position, and T arbitrary. But it 
is easy to see that by a re-definition of the space and time coordinates, we 
can make all 5x n vanish, keep g^y equal to zero, and make Sg^j equal to 

which in the limit of large a(t) approaches the time-independent function 
2(5jjC( x )) with a correction of order a~ 2 . In the non-linear case there is 



10 



no advantage to using something like Newtonian gauge (for instance, by 
choosing space coordinates so that g^ = e 2 ^[e 7 ]y where dijij = as well 
as 7jj = 0, as in [3]), and such a choice has the disadvantage of spoiling 
three-dimensional coordinate invariance. 

This analysis allows us to make a rough estimate of the expected cor- 
rections to the constancy of the metric correlations and to the vanishing 
of the correlation functions for matter perturbations following single-field 
inflation. At horizon exit the rate of change of the correlation functions of 
Qij is of order H, and from then until the end of inflation, the factor 1/a 2 
decreases by a factor roughly of order e~ 120 to e _140 [14], so at the end of 
inflation we expect the correlation functions of g^j to be changing at a rate 
of order e~ 120 H to e~ l ^°H. Similarly, the correlation functions of g^ at the 
end of inflation are of the same order as at horizon exit, so with the decrease 
in 1/a 2 we expect correlation functions for matter perturbations after infla- 
tion to less than the correlation functions of by a factor of order e -120 
to e~ 140 . This is the suppression factor that has to be overcome in order 
for physical processes during reheating like parametric amplification[12] to 
produce significant changes in observable correlation functions. 

III. Explicit Solutions for Multiscalar Theories 

The arguments of the previous section were rather abstract, so to see in 
more detail how they work out in practice, let us consider a more concrete 
but still fairly general model. To represent the matter fields, we suppose that 
there is a set of scalar fields Xm with a conventional kinematic Lagrangian 
and a completely arbitrary real potential V(x)- The unperturbed values of 
these fields are functions x n (t) of time that satisfy the field equations 

Xn + 3H Xn + ^^ = (5) 
where (in units with 8ttG = 1) 



ZH 2 = Wx n + V{x). (6) 



2 „ 



For the metric, we use the ADM parameterization[15] 
g 00 = -N 2 + g^WN* , g 0i = gijW = N> 

g 00 = -N- 2 , g i0 = N i /N 2 , g ij = {3) g ij - N*N j /N 2 , (7) 
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where ^g l i is the reciprocal of the 3 x 3-matrix g$j. It will be convenient 
also to write 

gij(x,t) =a 2 (t)g ij (x,t) , (8) 

where a(t) is the Robertson- Walker scale factor, satisfying a/a = H, with 
H given by Eq, (6), The quantities iV and ./V* are auxiliary fields, whose 
time derivatives do not appear in the Lagrangian. The Lagrangian for this 
theory is 

L=\jd 3 x y^Det^j-^ 4 ) - sTY.d^XnduXu ~ 2V(x) j 
= y Jd 3 x TV^j - a-tgVRij + C l 3 Ch - (C\) 2 
+ N- 2 {Xn ~ N'diXnf ~ a~ 2 p Yl diXndjXn ~ 2V(x)\ , 

n n ) 



(9) 



where <7 lJ (x,i) is the reciprocal of the matrix gijipc, t); g(x.,t) is the deter- 
minant of % (x, t) ; Rij (x, t) is the three-dimensional Ricci tensor (with the 
sign convention of [2]) for the metric ^(x, t); and C l j(x, t) is the extrinsic 
curvature of the surfaces of fixed time 

C l 3 = a- 2 ~g lk C kj , dj = ^ [2aa~ gij + a% - V t N, - VjN t ] , (10) 

where Vj is the three-dimensional covariant derivative calculated with the 
three-metric gij. For future use, we also note the well-known relations (for 
8ttG = 1): 

J2f n = -2H V(x) = 3H 2 + H. (11) 

n 

Models of this sort can be used both as fairly realistic theories of inflation, 
and also as surrogates for a theory of the matter and radiation after inflation. 
Because we are allowing any number of scalar fields, this model will in 
general have solutions in which neither the perturbations to matter fields 
nor the rate of change of g^ vanish at late times, so it is not trivial to see 
that there is also an adiabatic solution in which they do go to zero at late 
time, and that this is the solution that is excited if during inflation there is 
only one non-negligible scalar field. 

The gravitational field equations derived from this Lagrangian are 

V t (C^-5 l 3 C k k ) =±J2 d 3Xn(xn-N*d tX n) , (12) 

n 
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N 2 [ - oT^Rij - &pi + {C l if - 2V( X ) 

= J2(Xn ~ N l d lXn ? + N 2 a- 2 ~ 9 v £ d iXn djXn , (13) 



Rij - C k k C lj + 2C ik C k j +N~ 1 (- Cij + C\V 3 N k + C k 3 V l N k 
+ N k V k Cij + ViVjJv) = -a%F(x) - J] $Xn d jX n , (14) 

n 

and the scalar field equations are 



1 d ( \ d (y/gW 



dxi \ N 



dV(x) 

dXn 



(15) 



In line with the remarks of the previous section, we look for a solution in 
which <5% n (x, t) = x n (x, t) — x n (t) as well as and 5N = N — 1 are all of 
order a~ 2 (t) at late time. For convenience, in accordance with remarks at the 
end of the previous section, we also adopt a definition of space coordinates 
for which N l = 0. We can then write 



( " j - + Cj , 



(16) 



where Cj is> like and SN, a quantity whose leading term is of order a 2 : 



(17) 



Cj = ~ 2 9 lk [hj - lH5N~g k j\ + 0(l/a 4 ) . 
Then the gravitational field equations (12)-(14) become 

Vi (?j ~ S'^k) = dj J2 SXnXn + 0{a~ A ) , 

n 

AH5N = -a- 2 ~g^R l3 + AH^k ~ 2 ]T % n 5 Xn ~ 2 ]T ^^-S Xn + 0(a 

n n °Xn 

C 3 + mC j + HV 3 e k - HSNFj = a~ 2 g ik R kj + V 3 £ ^^-6 Xn + 0(0 . 



(18) 



(19) 



(20) 
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Using Eqs. (19) and (5), we can rewrite Eq. (20) in the form: 



+ 3HE^ = 1 \g ik R kj - \s l fg M R kl ] + 0( a - 4 ) , 



(21) 



where 



(22) 



CI 



Also, Eq. (18) now reads simply 



V i (H* i -^-H fc fc )=0(o- 4 ) 



(23) 



Eq. (21) has a solution: 



^•(x,t)= a ifc (x)^(x) 




(24) 



where T is any fixed time, Gij(x) is the value of gij (x, t) at that time, 1Zij(x) 
is the Ricci tensor calculated from the 3-metric ^ (x) , and B l j (x) is a time- 
independent function of co-moving position. It is convenient to choose T at 
around the end of inflation, where k/aH is smallest, in which case all terms 
in Eq.(24) are very small from soon after horizon exit to just before horizon 
re-entry. The leading term in (24) automatically satisfies Eq. (23) because 
of the Bianchi identity satisfied by IZij. 



To solve for the metric, we need to complete our choice of gauge. By 
using Eqs. (22), (17), (19), and (5), we have 



tij = 2g ik ~ k j + — n -~g i -E k k - -^~ gij ~ g kl R kl - ~ gij X + 0{a 4 ) , (25) 



where X = 0(1/ a?) arises from the perturbation to the scalar fields 



Under a shift t — > t+e(x, t) in the time coordinate, with e of order 1/a 2 (and 
a corresponding transformation x % — > x % + Q 1 ^ j dta~ 2 de/dxi to keep iVj = 
0), the perturbations to the scalar fields undergo the gauge transformation 
<5% n (x, t) — > <5% n (x, t) — e(x, £)x n (£) to order 1/a 2 . Hence, using Eq. (11), to 
this order 



2H 2 



H 




(26) 




(27) 
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so we can evidently choose e to make X = 0. This choice provides the great 
advantage that we can solve Eq. (25) for without first solving the field 
equations for the matter fields: 



9ij(*,t) = Qij{yL) + 2 



1 



7£y(x)- -^(x)g fci (x)^ fc/ (x) 



dt' 

t aW) Jt 



a(t") dt" 



1„ , ^ , x H 2 (t')dt 

+ -g ij {*)g k \*)n kl {yL) 

- ^(x)^(x)^(x)^ 



' ft' 



t H(t')a 3 (t') Jt 
1 H(t')dt' 
It a 2 {t')H{t') 



a(t") dt" 



V ; JV > Jt cfi{t>) JV ; V ' Jt a?{t')H{t') 
+ 0(a-\t)) , 



(28) 



where T is again some fixed time, conveniently chosen as the time at the 
end of inflation, and Qij{x) and Tl^j{x) are the values of gij and the asso- 
ciated Ricci tensor at that time. This confirms that while far outside the 
horizon, the time-dependent part of gij is of order a~ 2 . But in the radiation 
or matter-dominated era the second, third, and fourth terms in Eq. (28) in- 
crease like l/a 2 H 2 , which produces the breakdown in these approximations 
when physical wavelengths re-enter the horizon. 

It remains to consider the scalar fields. By using Eq. (5) again, we can 
put the field equation (15) in the form 



d 2 V(x) 

dXnOXm 

With X = 0, we now have 



&Xr, 



--g'tgij + 5N + 6H5N ) ^+25A^ n +0(l/a 4 ) 



SN = ~^~9 ij Rij + f + ^ £ *n *Xn + 0(l/a 4 ) . 



(29) 
(30) 



Using Eqs. (21), (25), and (30), we can put Eq. (29) in the form 



<%i + 3H5Xn ~ 7777 £ XnXm^X 



+ £ 



2H 

d 2 V{x) 

dXndXm 



H 

2H 2 



XnXr, 



2H 



= jf 2 {HX n ~ 2H Xn ) (j-z^j ~ HZ k ^ + 0(a" 4 ) , 



&Xm 
(31) 
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or, using Eq. (24) for S: 



1 



5Xn + 3H5Xn -W7tJ2 XnXm^X 



2H 

d 2 V(x) 



H 



+ ^ dxndxm \' 2J1 2 
= -k(Hx n -2Hxn 



XnXr, 



H 2 
+ 0(a~ 4 ) 



2H 

t 



(.XnXm ^XnXr 



$Xr. 



(32) 



An inhomogeneous differential equation of this form will have a solution 
in which a non-zero curvature scalar 1Z will generate perturbations of order 
1/a 2 in the various scalar fields, as anticipated in the previous section. The 
field equations also have isocurvature solutions in which 1Z = and there 
are small perturbations to the scalar field, not necessarily of order 1/a 2 , for 
which 



= SXn + 3H5Xn ~ XI XnXm$X 



+E 



d 2 V(x) 

dXndXrn 



- 3- 



m 

H s 
iff 2 



XnXm (-XnXm ^XnXm) 



$Xr, 



(33) 



and X = 0. To tell what solutions actually describe the metric and matter 
of the universe, we need a specific model of inflation, such as single field 
inflation, to which we now turn. 

IV. Single-field Inflation, and its Aftermath 

During single-field inflation there is by assumption only one non-zero 

■ 2 " ... 

Xn, say xi, so Eq. (11) gives H = — Xi/2 and H = — XiXi> an d we see that 
in this case the right-hand side of Eq. (32) vanishes. Thus during single- 
field inflation Eq. (32) is a homogeneous differential equation for 5xi, and 
therefore allows a solution Sxi = 0, which of course it must, since we can 
arrange that <5xi = by a choice of gauge consistent with the gauge choice 
X = used to derive Eq. (32). 

This shows that the non-linear field equations for single-field inflation 
have a solution in which 5xi = 0, and in which for late times <7y(x, t) is 
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attracted to a time- independent metric Qij(x), with corrections of order 
aT 2 and a~ 3 given by Eq. (28). We know by explicit calculation that in 
the linear approximation all solutions are in the basin of attraction for this 
asymptotic solution[16], but it is difficult to show that the relevant solution 
of the full non-linear equations is in this basin of attraction, and we shall 
simply assume that this is the case. 

Then at the end of inflation the transfer of energy from the inflation turns 
on other scalar fields, and the right-hand side of Eq. (32) becomes non-zero. 
As we have seen the general solution for the scalar field perturbations is a 
forced term of order a~ 2 , plus a solution of the homogenous equation (33). In 
general the solution of the homogeneous equation could be of any order in a, 
but by definition during single field inflation in our gauge all b\n and 5\ n are 
negligible, and with these initial conditions the solution of the homogeneous 
equation must be of order aT 2 to cancel the 0(a~ 2 ) terms in the solution of 
the inhomogeneous equation immediately after single field inflation. Thus 
as expected, for this solution all perturbations to the matter fields become 
of order 1/a 2 outside the horizon, and we have a pure adiabatic solution, 
with negligible corrections. 

V. Tree- Approximation Correlation Functions 

If the results we have obtained so far really applied to the metric and 
matter perturbations in the Heisenberg picture, we could conclude that with 
a suitable definition of coordinates, all correlation functions involving only 
Qij (or functions of g^j) become time- independent outside the horizon, and 
that all correlation functions involving perturbations to goo, 9oi, an d mat- 
ter variables become negligible outside the horizon. But as mentioned in 
the Introduction, the presence of quantum fluctuations of arbitrarily small 
wave lengths invalidates the expansions in powers of 1/a as applied to the 
Heisenberg picture interacting fields. To avoid this problem we must limit 
our consideration to tree graphs for correlation functions, on the assumption 
that the contributions of graphs with loops are much smaller. We can as 
usual apply the results of Sections II - IV to the Heisenberg picture quantum 
fields, but calculate correlation functions only to lowest order in interactions 
to avoid loop graphs, hoping that this is a good approximation. Here we 
want to consider an alternative approach, in which one explicitly considers 
only tree graphs. 

In the appendix we review the general tree theorem of [11], which shows 
how to calculate the sum of tree graphs for correlation functions by a so- 
lution of the classical field equations, subject to certain constraints. To 
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illustrate the use of this theorem, in this section we will apply this theo- 
rem to the correlation functions of the reduced metric g^ = gij /a 2 during 
single-field inflation, adopting space and time coordinates for which there is 
no perturbation to the inflaton field, and for which gio = 0. 

The generating function W[J, t\] for correlation functions of g^ at a time 
t\ is defined by Eq. (A.l), which for this case takes the form 



exp{VF[J,ti]} = ^0,in exp [ J d 3 x g? (x,h)J ij (*) 



0,in> , (34) 



where gfj (x, t) is the Heisenberg-picture quantum mechanical operator cor- 
responding to (jij(x,t). Correlation functions for gij(x., ii) are calculated 
according to Eq. (A. 2), which here reads 



(0, in|^(x,*i)^(y, ti)--- |0, in) 



0J«(x) d,J kl (y) 



■ exp 



{w[J,h}} 



(35) 

We want to evaluate W[J, ti] for late times during inflation, at which the 
Robertson- Walker scale factor a(t\) becomes exponentially large, from which 
we can calculate the late-time expectation value of products of the operators 
gij at various space coordinates or wave numbers. 

As described in the appendix, to calculate W in the tree approxima- 
tion we construct complex c-number metric fields gij(x,t) together with a 
complex auxiliary field iV(x, t), satisfying the constraints: 

(A) The fields satisfy the Euler-Lagrange equations. In the our case, they 
are Eqs. (12)-(14) with no scalar field perturbations. 

(B) The fields g^ satisfy constraints at time t\\ 

(36) 

-J^'(x) . (37) 




(C) gij satisfies a positive frequency constraint at time t — > — oo, that it 
behaves as a superposition of terms proportional to exp(— iut) , with uj 
various positive frequencies. 

These constraints give the functions gij(x, t) an implicit dependence on both 
the current J and on the time t\ at which correlations are to be measured. 
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With the functions «7ij(x, t) and 7V(x, t) constructed in this way, the gener- 
ating function is given by Eq. (A. 6), which here reads 

W[J,h] tree = f 1 lmL[g{t)j(t),t]dt+ f d 3 x J^'(x)^ (x, h) . (38) 

We showed in Section III that the non-linear field equations have a solu- 
tion for gij(x, t) that for late times is attracted to a time-independent metric 
Qij (x). But as remarked in the appendix, this is not enough to conclude that 
the correlation functions for gij(x,t) become time-independent at late time. 
We must also show that the constraints (36) and (37) do not give 5ij(x, t) 
any dependence on the time t\ at which the constraints are imposed, pro- 
vided a(t\) is sufficiently large, and we must consider the convergence of the 
time integral in Eq. (38) as a(t\) at the upper limit t\ becomes large. 

For large a(t\), the constraint (36) simply provides the ti-independent 
condition that the leading term Qij(x) in Eq. (28) is real for all x. It follows 
then that the Ricci tensor 7£jj(x) for the metric Qij(x) is real, so the terms 
in (28) of order I /a 2 are also real. The leading terms in Im^(x, t\) are 
then of order 1/a 3 : 

ft Afl ft ff2(t')f]t> 

yjK ' ' K ' JK ' Jt a 3 (t') ' kK ' Jt a?(t')H{t>) 

+ 0(a-\t)). (39) 

The functional derivative appearing in the constraint (37) is 



iMMM = « W f^ (x, t )(-2i f ( t )y t+ ^, t )-^',(*, t )) 

(40) 

The metric is constrained by Eq. (36) to be real at t = t\, so the term 
—2H5' l j in parentheses makes a contribution to this functional derivative 
that is also real at t = t±, but the tensor £ l j in the other two terms has 
an imaginary part given by the 0(a -3 ) term in Eq. (24) (with 3 l j replaced 
with which in the absence of scalar field perturbations is the same): 

Imf j-(x,t) = a _3 (t)Imfi i 3 -(x) + o(a" 4 (t)) , (41) 

so 



(42) 
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Thus the constraint (37) does become independent of t\ for large t±. 

This is not just a happy accident. We can understand the asymptotic 
constancy of the left-hand side of Eq. (37) by recalling the Euler-Lagrange 
equations 

d 5L[g{t)M)\ _ SL[g(t),9(t)} 
dt <5g i3 -(x,i) <%(x,t) 

The imaginary part of the right-hand side decreases as I /a 2 , so the left-hand 
side of Eq. (37) becomes constant for large a. 

This leaves the question of the convergence of the integral over time in 
Eq. (38) for large a(t\). Let's first consider the terms in the gravitational 
part of the Lagrangian (9) that contain either or 1 space or time derivative. 
Since N is fixed by the condition that the Lagrangian be stationary in N, to 
first order in 5N we can set N = 1. It is then straightforward to calculate 
that the terms in the gravitational part of L of zeroth or first order in 
derivatives add up to 

Li[g,~g, t] = y J d 3 x^ [ - 12H 2 -AH- 2#2% - SHgVViNj] . (43) 

The final term in square brackets integrates to zero (and in any case vanishes 
for the choice we have made of spatial coordinates), leaving us (as already 
noted in [3]) with a total time derivative 

Lx [g, ~g, t] = ~2j t (a 3 H J d 3 x^j (44) 

As remarked in the appendix, such a total time derivative in the Lagrangian 
has no effect on the correlation functions. This leaves the terms in the 
gravitational part of L that are of second order in £ or that involve the 
space curvature. According to Eq. (39) the imaginary part of N \fgg i:> Rij 
is of order a(i)~ 3 at late time, which cancels the over-all factor a 3 in the 
Lagrangian, so this term makes a contribution of order a~ 2 . According to 
Eq. (41), the imaginary part of any second-order function of the £ is of 
order a~ 2 x a -3 , so again such terms make contributions to Im(L — L\) 
that at late times are of order a~ 2 . The time integral in Eq. (38) therefore 
converges to a finite limit for large a{t\) exponentially fast, as J* 1 a{t)~ 2 dt. 
This concludes the proof that in single field inflation the generating function 
W[J, ti] converges to a ti-independent function for large t±, and therefore 
so do the correlation functions of g^. 
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This demonstration, that the correlation functions for the metric con- 
verge to t\ -independent functions for large t±, does not imply that these 
limits are uniform in the coordinates appearing as arguments of the metric 
components. In fact, if we set coordinates equal, the correlation functions 
do not converge to finite limits. For instance, if to avoid ultraviolet diver- 
gences we define ( r (t) as the average of the curvature perturbation £(x) over 
a very small co-moving volume r 3 whose physical radius a(t)r over the times 
of interest remains much less than the horizon size 1/H(t), then it can be 
shown that in slow roll inflation the tree-approximation vacuum expecta- 
tion value of (r(t) increases like a(t) ras_1 for n$ > 1 (where ns is the usual 
scalar mode slope parameter) and like lna(t) for ng = 1, though it does 
approach a constant for rag < 1. Because of the way that a(t) and the co- 
moving coordinates x 1 enter in the flat-space Robertson- Walker metric, they 
have no physical significance in themselves; it is only a(t) times differences 
of co-moving coordinates that have a significance, as physical separations. 
Thus we expect the metric correlation functions to approach constant limits 
only when all such physical separations become large compared with the 
horizon size 1/H. Of course, in practice we are chiefly interested in the 
Fourier transforms of the correlation functions, in which case the physical 
wave numbers are the co- moving wave numbers divided by a(t), and we ex- 
pect these Fourier transforms to approach finite limits only when all physical 
wave numbers become much less than H(t). 

I am grateful for helpful conversations with Raphael Flauger, Eiichiro 
Komatsu, David Lyth, Juan Maldacena, Misao Sasaki, and Richard Woodard 
This material is based upon work supported by the National Science Foun- 
dation under Grant No. PHY-0455649 and with support from The Robert 
A. Welch Foundation, Grant No. F-0014 

Appendix: The Tree Theorem 

In this appendix we will review the general tree theorem of [11], in a 
somewhat simplified notation, and add a remark that is needed in Section 
V. We consider a general Lagrangian system, with Hermitian Heisenberg- 
picture canonical operators q^it), and Lagrangian (not Lagrangian density) 
L[q H (t),q H (t),t], possibly with an intrinsic time dependence. In field theo- 
ries the index a incorporates a space coordinate x as well as discrete indices 
labeling the various field components; a sum over a includes an integral over 
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x as well as sums over discrete indices; and derivatives with respect to q a (t) 
are interpreted as functional derivatives. We wish to calculate the generat- 
ing function W[J, t±], a real function of a set of real c-number currents J a , 
which is defined by 



exp{w[J,t 1 }} = /o,in exp [j>*(ti)J a 



0,in 



(A.l) 



where |0,in) is a state defined to look like the vacuum state at an early 
time, which in this paper we take as t = — oo. From W we can calculate 
expectation values in this state of products of any number n of q H s at the 
time t\. 



(0, inlgf (tx) gf (*!)••• |0, in) = 



gn 



dJ a 9J b 



■ exp 



{W[j,h]} 



j=o 



(A.2) 



To calculate W[J, t±] in the tree approximation, we construct complex c- 
number functions q a (t), subject to three conditions: 

(A) The q a {t) satisfy the Euler-Lagrange equations 



8 dL[q(t),q(t),t] _ dL[q(t),q(t),t] 



dt dq a (t) 



dq a (t) 



(A.3) 



(In extending the Lagrangian to complex variables, we take it as a real 
function, in the sense that L*[q(t),q(t),t] = L[q*(t),q*(t),t\.) 

(B) The q a (t) satisfy constraints at time t\: 



Img a (ii) = , 



Im 



9L[g(ti),g(ti),ti] 
dq a (ti) 



= ~Ja 



(A.4) 
(A.5) 



(C) The q a (t) also satisfy a positive frequency constraint at time t — > 
— oo, that they behave as superpositions of terms with time-dependence 
exp(— iut), with u various positive frequencies. 

(In [11] the functions q a (t) were denoted qLa(t); we are here taking advantage 
of the fact that for real currents, the other functions qR a (t) introduced in [11] 
are just q*L a {t)-) The constraint (B) gives the q a (t) an implicit dependence 
on t\ as well as on the J a . With q a (t) calculated subject to these three 
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constraints, the contribution of connected tree graphs to the generating 
function is given by 

W[J,h] tiee = f 1 Im L[q(t),q{t),t}dt + Y,Ja<la(h) . (A.6) 

J-oo a 

We are concerned in this paper with the limit of W for large t\ (or more 
precisely, for large a(t\)). 

Prom the foregoing, we can see that, in order to conclude that the gen- 
erating function becomes independent of t\ when t\ is sufficiently large, it 
is not enough to show that the quantities q a (ti) approach finite limits for 
large t\. We must also show that the integral in Eq. (A.6) converges in 
this limit. (There is no problem with the convergence at very early times, 
where the integrand oscillates increasingly rapidly.) Further, because the 
constraints (A. 4) and (A. 5) are applied at time t\, we must show that the 
quantities lm{d L[q(t\), q{ti), ti}/dq a {t\)} as well as Imq a (ti) approach finite 
t\ -independent limits for large t\. 

In order to evaluate the late time behavior of the correlation function in 
Section V, we need to supplement this general review with a remark about 
the effect of adding to the Lagrangian a derivative term: 

AL(t) = j t F[q(t),t] , (A.7) 

with F[q(t),t] an arbitrary function of t and of the q a (t), which is real in 
the same sense as L; that is, F*[q(t),t] = F[q*(t),i\. It is familiar that 
such derivative terms do not matter in calculating the S-matrix, because 
there the Lagrangian enters in integrals over all time, but in calculating the 
generating function here we need to integrate the Lagrangian only up to time 
ti, and the Lagrangian also enters in the constraint (A. 5). Nevertheless, we 
can easily see that in calculating the generating function, as in S-matrix 
calculations, the change (A.7) has no effect. First, adding a derivative term 
(A.7) obviously has no effect on the Euler-Lagrange equations (A. 3). The 
only other place where the Lagrangian enters in constructing the functions 
q a (t) is in the constraint (A. 5), but adding the derivative term (A.7) changes 
the left-hand side of Eq. (A. 5) by 

A Im MfaM>».>.'.] _ Im ^MlMil , ( A .8, 

dq a (h) dq a (ti) 

and this vanishes because the constraint (A. 4) requires that q a (ti) be real. 
Hence the change (A.7) has no effect on the functions q a (t). The only effect 
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on the generating function (A. 6) is then to change it by an amount 

rti 

AW[J,h] tiee = / 1mAL[q(t),q(t),t]dt = 1mF[q(t 1 ),t 1 ] (A.9) 

J — oo 

and this vanishes because again the constraint (A. 4) requires that q a (ti) be 
real. 
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